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We argue that atmospheric cascades can be regarded as example of the self-organized criticality and studied by 
using Levy flights and nonextensive approach. It allows us to understand the scale-invariant energy fluctuations 
inside cascades in a natural way. 


1. INTRODUCTION 

It is well known that energy spectra of par¬ 
ticles from atmospheric family events which are 
observed by the emulsion chambers at mountain 
altitudes are essentially following power-like de¬ 
pendence. On the other hand the occurence of 
power-law distributions is a very common fea¬ 
ture in nature and it is usually connected with 
such notions as criticality, fractals and chaotic 
dynamics and studied using generalized (nonex¬ 
tensive) maximum entropy formalism character¬ 
ized by a nonextensivity parameter q. Such for¬ 
malism leads in a natural way to power-laws in 
the frame of equilibrium processes [11]. For non¬ 
equilibrium phenomena the sources of power-law 
distributions are self-organized criticality Q and 
stochastic multiplicative processes |. In the 
former nonequilibrium systems are continuously 
driven by their own internal dynamic to a critical 
state with power-laws omnipresent. In the later 
power-law is generated by the presence of under¬ 
lying replication of events. All these approaches 
can be unified in terms of generalized nonexten¬ 
sive statistics mentioned above (widely known as 
Tsallis statistics) jli]. 

We shall look therefore at the development of 
cascades from this point of view. To be more 
specific, we shall use the notion of Levy flights 
as representatives of power-laws emerging from 
generalized statistics [Q. 


2. MONTE-CARLO SIMULATION OF 
MULIPLICATIVE PROCESSES 



Fig.l Distributions P(x) for multiplicative 
process, as given by eq. (jl[) (histograms), for 
different generation numbers N = 2,6 and 
8 are fitted by Levy distribution (eq. (|j)) 
(solid lines) with parameters a = 4.81,2.05 
and 2.015, respectively. 

Let us start with numerical (Monte-Carlo) con¬ 
sideration of the following model of multiplicative 
process: 


x 0 = 1 , 







2 


where £ is random number taken from the expo¬ 
nential distribution (to account for the essentially 
exponential scaled energy dependence of the sin¬ 
gle elementary interaction): 


/(£) = h exp (-£) • <2) 

As was shown in [Q, every exponential distribu¬ 
tion with a fluctuating parameter results in the 
Levy distribution (cf. Fig. 1), which in our case 
has the following form 


P{x N ) 


1 a - 1 
(xn) a — 2 


1 x N 
a — 2 (xjv) 


.(3) 



where c is generation independent parameter 
which should be fitted (in our case c = 0.55). 

If we use the canonical form of Levy distribu¬ 
tion as emerging from Tsallis statistics 


P(x) = 


2 -q 

X 




( 5 ) 


where mean (x) is given by the parameter \ = 
(3—2 q)(x) and q is mentioned before nonextensive 
parameter such that q = 1 + 1/a j| then 


Q = 


3 c N ~ l 
2 2 


( 6 ) 


Notice that for IV —> 1 parameter q —> 1 (or, 
respectively, a —> oo), i.e., we are recovering 
the initial exponential distribution. In the limit 
N —> oo parameter q approaches value q = 3/2, 
which is limiting value available for q emerg¬ 
ing from the normalization condition imposed 
on the probability distribution P(x at). In this 
limit a —> 2. It is interesting to note that for 
[xn/(xn)\-[^/{ ot — 2)] >> 1, i.e., for a; sufficiently 
large one gets power-like behaviour of P(x at): 


P(M) “ (<Sy) ■ (7) 

Actually such situation is reached reasonably fast, 
because (xjv) = and a tends to its limiting 
value a = 2 rather quickly with increasing num¬ 
ber of generations N. In practice the equilib¬ 
rium distribution P(xn ) oc xjj 2 is reached (for 
xn » (xn)) already for N > 6. 


Fig.2 Dependence of the parameter a of the 
Levy distribution (|3[) on the generation num¬ 
ber N. Solid line represents fit by the for¬ 
mula (|l|) with c = 0.55. Asymptotically 
a = 2 (dotted line). 

In our present case parameter a depends on the 
number of generations N and on the parameter 
£o which defines exponential distribution of the 
variable £ (in the limit of large N dependence on 
£o effectively vanishes). As is illustrated in Fig. 
2, dependence a(N) can be described by simple 
formula 


3. COMPARISON WITH EXPERIMEN¬ 
TAL DATA 

There exists a number of experimental data 
from the emulsion chambers exposed at mountain 
altitudes, which are relevant for our approach and 
which we shall compare with now ||,[| . Although, 
as we have already mentioned, the energy spectra 
of particles from atmospheric family events they 
represent are roughly expressed by the power-law 
distributions, so far there was no analysis show¬ 
ing wheather it is true and to what extend. Our 
work is the first attempt in this direction as we 
show that they all can be described by the Levy 
type spectra. 


1 - c N ~ l 
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Fig.3 Integral energy spectra for gamma 
quanta (circles) and hadrons (squares) in 
families registered at Mt. Kambala are fit¬ 
ted by Levy distributions (solid lines). 

Integral spectra TV(> E) from Kambala experi¬ 
ment (at 520 g/crn 2 ) [|| are shown in Fig. 3 where 
they are fitted using parameters a — 1 = 1.8 (for 
gamma families) and a — 1 = 2.0 (for hadronic 
families). Notice that hadronic component is lit¬ 
tle ’’younger” (higher a means smaller TV, cf. Fig. 
2) than the electromagnetic component. 

Families observed deeper in the atmosphere 
(what means higher TV) have smaller a, as ex¬ 
pected, cf. Fig. 4, where N(> E) from Pamir ex¬ 
periment (at 600 g/cm 2 ) j(ij] are presented. Here 
the corresponding parameters a — 1 are equal to 
1.8 and 1.7 for gamma and hadronic components, 
respectively. 

It is evident from Figs. 3 and 4 that the ob¬ 
served gamma-hadron families are already reach¬ 
ing their quasi-equilibrium states. It should be 
noticed that most of the families registered come 
not from a single nuclear interaction but they usu¬ 
ally contain particles which are decascadents of 
several (~ 7 at E E > 100 TeV) genetically con¬ 
nected nuclear interactions [[ij. Our analysis in¬ 
dicates that the average number of cascade gen¬ 
erations leading to the observed distributions is 
about TV ~ 3. 
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Fig.4 Integral energy spectra for gamma 
quanta (circles) and hadrons (squares) in 
families registered at Mt. Pamir are fitted 
by Levy distributions (solid lines). 


4. SUMMARY AND DISCUSION 

We have provided description of cascade pro¬ 
cesses encountered in cosmic ray emulsion exper¬ 
iments at high altitudes, which is based on the no¬ 
tion of Levy distrubutions (which, in turn, orig¬ 
inate from the nonextensive statistics described 
by parameter q jl]|4|). In this way we have found 
that the observed distributions are to a high ac¬ 
curacy power-like, but at the same time we are 
able to account for the small deviations from the 
exact power-like behaviour. They are, in our ap¬ 
proach, directly connected with the finite number 
TV of the cascade generations in the way provided 
by formulas (|j) or (j(i). 

The results presented before should be con¬ 
fronted with the result of pure nuclear origin, 
which is the case of the so called ’’halo” events 
with strong concentration of particles in the cen¬ 
tral region. Fig. 5 shows such a case expos¬ 
ing integral energy distribution of shower cores 
(recorded at distance R < 0.65 cm from the 
energy-weighted center) obtained in a ’’halo” 
event P06 registered at Chacaltaya ||. In this 
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case we observe almost exponential distribution 
(a — 1 = 16.0 in this case, what translates to 
q = 1.06) with no influence of multiplicative pro¬ 
cesses. 



Fig.5 Integral energy spectra for shower cores 
from halo event (circles) registered at Cha- 
caltaya is fitted by Levy distribution (solid 
lines). 

Finally, it should be noticed that there exists a 
number of special cases more-or-less directly con¬ 
nected with our results: 

(а) Energy transport equation 0 for cascade pro¬ 
cesses gives 

N(E, t) = N 0 E -( s+1 > exp • f) (8) 

v s +1 t 0 J 

as number of particles with energy E at depth t. 
Notice that for the age parameter s —> 1 one ap¬ 
proaches equilibrium spectrum N(E) = No E~ 2 
with the slope a = 2. 

(б) One observes strong intermittent behaviour in 
the families (jlo), which results from the fluctua¬ 
tions in the atmospheric cascades alone and is not 
sensitive to details of the elementary interactions, 
(c) One observes multifractal behaviour in the 
cascades with fractal dimension at each cas¬ 
cade stage N where P(x) oc x Dn Jtl). For 
each successive emission the distribution should 


be more inhomoheneous, i.e., Djy > Dn + 1 . 

(d) Self-organized criticality to which all cas¬ 
cades considered here bear close resemblance be¬ 
cause of the Zipf’s law iQ , which predicts power- 
law behaviour of the type seen in eq.(Jj]) (with 
a = 2 ). 
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